Stochastic Electron Acceleration in Shell-Type Supernova 

Remnants 

Siming Liu^, Zhong-Hui Fan^'^, Christopher L. Fryer^''^, Jian-Min Wang^'^, and Hui Li 



ABSTRACT 

We study the stochastic electron acceleration by fast mode waves in the tur- 
bulent downstream of weakly magnetized collisionless astrophysical shocks. The 
acceleration is most efficient in a dissipative layer, and the model characteristics 
are determined by the shock speed, density, magnetic field, and turbulence decay 
length. The model explains observations of shell-type supernova remnants RX 
J1713. 7-3946 and J0852. 0-4622 and can be tested by observations in hard X-rays 
with the HXMT and NuSTAR or 7-rays with the GLAST. 



Subject headings: acceleration of particles 
— turbulence 



MHD — plasmas — shock waves 



Introduction 



Recent observations of high-energy emission from a few shell-type supernova remnants 
(STSNRs) advance ou r understanding of the underlying physical processes significantly 
(lEnomoto et al. 1120021). These observations not only give high quality emission morphology 



and spectra 



Hiraga et al. 



Jchivama et al. 



J cnivi 
2OO5I : 



2003: 



Takahashi et al. 



Aharonian et al. 2004; Cassam-Chenai' et al. 2004 



2008al ). which lead to tight upper limits on thermal 



X-ray s, but also discover X-ray variability on a timescale of a few months (jUchiyama et al. 



20071 ). which is intimately connected to processes near the shock front (SF). The radio to 
X-ray emissions are produced by relativistic electrons through the synchrotron processes. 
The nature of the TeV emission is still a matter of debate. The major challenges to the 
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hadronic scenario include: 1) the lack of correlation between the TeV emission and molec- 
ular cloud distribution; 2) the requirement of very high mean target gas density, supernova 
explosion energy, and/or proton acceleration efficiency, which implies efficient magnetic field 
amplification in the cont ext of the diffusive shock acceleration and strong suppressio n of 



the electron acceleration (ILucek fc Bell Il2000l : iBerezhko fc Volk 



the good co rrelation between the X-ray and TeV emission (jPlaga I l2008l : lAharonian et al. 



2006 



Butt et al. 



2008); 3) 



20061 . l2007al ). which suggests a leptonic scenario. Because secondary leptons from hadronic 
processes won't produce as much radiation power as 7-rays from neutral pion (vr") decays 
( lAharonian et al. Il2007bl ). the fact that the X-ray luminosity is slightly higher than the TeV 
luminosity requires that most leptons be accelerated from the background plasma. 

Compared to the hadronic scenario, the leptonic scenario has much fewer parameters 
and is well-constrained by the observed radio. X-ray, and TeV spectra. In the simplest 
one-zone model, the TeV emission is produced through the inverse Compton (IC) scattering 
of the Galac tic interstellar radia tion by the same relativistic electrons producing the radio 
and X-rays (jPorter et al. II2006I ). The TeV spectrum can be used to infer the high-energy 
electron distribution. One then needs to fit both the synchrotron spectral shape and flux 
level by adjusting the magnetic field alone. TeV observations always se e spectral soften 



ing wi th the increase of the photon energy, suggesting a spectral cutoff (lAharonian et al 
200730). T here are also indications t hat the synchrotron sp ectrum cuts off in the hard 



X-ray band ( Takahashi et al. 2008al jbl; lUchiyama et al. " 2007). A single power-law model 
with an exponential cutoff for the electron distribution has difficulties in reproducing the 
broad TeV spectrum. Models with more gradual cutoffs lead to much better fits. Relativis- 
tic particle distributions with gradual high-energy cutoffs are a natural consequence of the 



stochastic acceleration (SA) by plasma waves in magnetized t urbulence (jPark fc Petrosion 



I995I : iBecker. Le fc Dermer I l2006l : IStawarz fc Petrosian 1 12008| ) . The particle distribution is 
determined by the interplay among the acceleration, cooling, and escape processes with the 
cutoff shape determined by the energy dependence of the ratio of the acceleration to cooling 
or escape rate. Given the low magnetic fields derived from leptonic models, the radiative 
cooling time is much longer than the remnant lifetime. The gradual high-energy cutoff has 
to be caused by the balance between the acceleration and escape processes. 

We consider the evolution of weakly magnetized turbulence in the downstream of strong 
non-relativistic shocks and study the SA by fast n iode waves. The large scale turbu lence 
cascades following the Kolmogorov phenomenology (IBoldyrev II2002I : iJiang et al. II2008I ) until 
it reaches the scale, where the eddy turnover time becomes comparable to the period of fast 
mode waves. The Iroshnikov-Kraichnan phenomenology prevails on even smaller scales. The 
coUisionless damping of plasma waves sets in at the coherent length of the magnetic field, 
where the eddy turnover time is comparable to the period of Alfven waves. Fast mode 
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waves propagating nearly parallel to the mean magnetic field can survive the transit-time 
damping (TTD) by the thermal background particles, and accelerate some electrons through 
cyclotron resonances to a power-law high-energy distribution, which cuts off at the energy, 
where the particle gyro-radius reaches the coherent length of the magnetic field (jhiu et al 



20061 ). Acceleration of higher energy electrons by the nearly isotropic fast mode waves on 



larger scales leads to a gradual cutoff. SA of electrons can naturally account for emissions 
from the STSNRs with a weak magnetic field. Future observations in the hard X-ray and 
MeV bands may test the model. The particle acceleration model is described in § [2] and 
applied to SNRs RX J1713-3946 and J0852.0-4622 to derive quantitative results in §[1 §11 
gives the conclusions. 



2. Turbulence Cascade, Wave Damping, and Stochastic Particle Acceleration 

The coUisionless nature of weakly magnetized shocks implies difficulties in converting 
the free energy in the system into heat and energetic particles instantaneously at a nar- 
row SF. Instead, the high magnetic Reynolds number implies strong turbulence. Previous 
studies of diffusive shocks focus on its roles in making high-energy particles crossing the SF 
repeatedly. The turbulence can also carry a significant fraction of the released free energy 
and decay gradually. The turbulence genera tion mechanis m is not well understood. It is 



usually attributed to the Weibel instability (IWeibel Ill959l ) of particle streaming and may 



also depend on the properties of the upstream plasma. We here assume that the turbulence 
is isotropic and generated at a scale L with an eddy speed u. The free energy dissipation rate 
is then given hy Q = Cipu^/L, where Ci ~ 1 is dimensionless. For strong non-relativistic 
shocks with the shock frame upstream speed U much higher than the speed of the upstream 
fast mode waves Vp = {v\ + 5f|/3)^/^, mass, momentum, and energy conservations across 
the SF require U"^ = 5f| + 5u^ + 2v\ + f/^/16 , where the Alfven speed va = (-B^/47rp)^/^, 
the isothermal sound speed vs = {P/pY^'^, and B, p, P are the magnetic field, mass density, 
pressure, respectively. The shock structure can be complicated due to the presence of strong 
turbulence, and the speeds vs, va, and u should be considered as averaged quantities over 
~ L, which also corresponds to the energy dissipation scale. The excitation of plasma waves 
should be very efficient in the supersonic phase with u > vp, and the waves may prevail in 
the downstream accelerating some particles in the background plasma to very high energies. 

The eddy turnover speed and time are given respectively by vl^^ik) = 4:7iW{k)k^ and 
Teddik) = 2tx /Cikvedd, where W{k) = {47r)-\2'KQ/Cipf/^k-^^/\ k = 27r/Z, and I are the 
isotropic turbulence power spectrum, wave number, and eddy size, respectively. At the 
turbulence generation scale km = ^rr/L, Vedd = u, and the total turbulence energy is given 
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by J W{k)4iTk^dk = (3/2)m^. The MHD wave period is given by TF{k) = 2Ti/vFk. Then the 
transition scale from the Kolmogorov to Kraichnan phenomenology kt occurs at Tp{kt) = 
Teddikt), which gives kt = {Ciu/vF)^km ■ For k > kt, the turbulence spectrum in the inertial 
range is given by 

W{k) = {A7r)-\vF/Cry/\'/'kl/'k-'/' . (1) 

The coUisionless damping starts at the coherent length of the magnetic field Id = 27r//c^, 
where the period of Alfven waves 2iT/kvA is comparable to the eddy turnover time, k^ = 
{Cfu^Vp/v\)km . 

Both va and vs may be lower than u near the SF. li u va ^ vs at the SF, fast, slow, 
and Alfven waves can all be excited, and the consequent plasma heating quickly leads to 
i^s ^ Va- Fast mode waves with the highest phase speed should dominate the SA processes. 
We assume that they dominate the energy dissipation as the turbulent flow moves away from 
the SF and vp evolves from va to u at the same time. For a fully ionized hydrogen plasma 
with isotropic particle distributions, whic h are reasonable in the abse r ice of strong large scale 



magnetic fields, the TTD rate is given by (IStix Ill962l : iQuataert Ill998l : iPetrosian et al. II2006I ) 



AT{9,k) 



{27rkBY/^ksmH 
2(me + m„) cos 6* 



(Teme)^''^exp 



mfUJ 



2kBT,kl 



+ (Tpmp)^''^exp 



rripUj 



where Tg, Tp, mg, m^, 6*, a;, and /c. 



(2) 

p, ""p, ^5^5 — COS 6* are the electron and proton temperatures, masses, 

angle between the wave propagation direction and mean magnetic field, wave frequency, 
and parallel component of the wave vector, respectively. The first and second terms in the 
brackets on the right hand side correspond to damping by electrons and protons, respectively. 
For weakly magnetized plasma with va < vs, proton heating always dominates the TTD 



for io'^/kn 



2k-QTp/mp (Quataert 1998). If va does not change dramatically in the 



downstream, the continuous heating of background particles through the TTD processes 
makes Tp — > {mp/me)Te since the heating rates are proportional to (mTY^'^, where m and 
T represent the mass and temperature of the particles, respectively. Parallel propagating 
waves with sin ^ = are not subject to the TTD. Obliquely propagating waves are damped 
efficiently by the background particles. Although the damping rates for waves propagating 
nearly perpendicular to the magnetic field w ith cos ^ ~ are also lo w, these waves are subject 
to damping by magnetic field wandering (IPetrosian et al. 1 120061 ) . The turbulence power 



spectrum cuts off sharply when the dampin g rate beconies com parable to the turbulence 
cascade rate F = T'^^/^Tp^ + t~2j) — t~^^tf (j Jiang et al. II2008I ). One can define a critical 



propagation angle Odk), where At{Oc, k) = T{k). Then for k^Tp = C2pTnpV? with C2p < 3/5 
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for vf < u, equations ([T]) and ([2]) give 



— )± 1^ exp I - 

2V27r3/2C'V2yt;^^l/2 COS V SCspM^ COS^ ^, 



sin^ 6c 



where the electron pressure and damping have been ignored. 



For va ^ we have 3f|./5 



[?>C2pllQ{C2p + The accelera- 



tion and scattering times of relativistic electrons by these nearly parallel propagating waves 
are energy independent and given, respectiv ely, by Tsc = 2Tx/ckd, and Tac = {'ic^/v%)Tsc 



( iLiu et al. I l2006l : iBlandford &: Eichler 1119871 ). where c is the speed of light. The electron 
escape time from the large scale eddies is given by t^sc = L"^ I^c^'Tsci and the spectral index 
of the accelerated electrons in the steady state is given by p = (9/4 + Tad'Tesc)^^'^ ~ 1/2 = 

1 /2 

(9/4 + Vl^v\IClu''v%Y'' - 1/2 . The maximum energy that electrons can reach is given by 
7cmec2 = 2'KqB/kd = 7.2Q[{p + 0.5)^ - 2.25Y/^[C2p/{C2p + 1)]^/^ {L/W^^cm){B /WnG)Uo 
TeV, where Uq = U /0.015c, q is the elementary charge units. The ratio of the dissipated 
energy carried by non-thermal particles to that of the thermal particles should be greater 
than T] = 9l{kd)/2 = e^^^v\/2T{{2T:C2pY^'^uvF, where e = 2.72, since the isotropic turbulence 
with k < kd can also accelerate particles with the Lorentz factor 7 > 7c. 

Electrons with even higher energy interact with the nearly isotropic fast mode waves 
at A; < kd- The corresponding acceleration and escape times are proportional to 7^/^ and 



7 



-1/2 



respe ctively, which leads to a steady-state electron distrib ution / oc exp— (7/7, 



\l/2 



for 7 ^ 7c ( jPark fc Petrosion Ill995l : iBecker. Le fc Dermer 1120061 ). To obtain the complete 
electron distribution from the thermal energy to the high-energy cutoff self-consistently, one 
needs to obtain the acceleration and scattering rates by these nearly parallel propagating 
(with k > kd) and nearly isotropic (with k < kd) fast mode waves and solve the particle 
kinetic equation numerically. These are beyond the scope of this paper. Instead, the above 
discussions indicate that the particle distribution may be approximated reasonably well by 
/ oc 7~^exp — (7/72)^^ with /3 = 1/2. To have significant particle acceleration, the acceler- 
ation time needs to be shorter than the turbulence decay time: Tac < Teddikm) = L/Ciu, 
which also validates the usage of the steady-state solution of the particle kinetic equation 
and implies Ci < {20C2p)^/y3[{p + 0.5)^ - 2.25]^/\ 



3. Results 

Several STSNRs have been observed extensively in the radio. X-ray, and TeV bands. X- 
ray observations with Chandra, XMM-Newton, and Suzaku, and TeV observations with HESS 
have made several surprising discoveries that challenge the classical diffusive shock particle 
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accel eration model. The SNR RX Jl 713. 7-3946 is about t = 1600 years old (IWang et al 



19971 ) with a radius of ~ 10 pc and a distance of D ~ 1 kpc. By fitting its broad- 



band spectrum with the above electron distribution and background photon field given by 



Porter et al. I (120061 ) through the synchrotron and IC processes, we find that p = 1.85, 



B = 12.0 /iG, 'yr.me.c = 3.68 TeV, and the total energy of electrons with 7 > 1800 (iLiu et al. 



20061 : iFan et al. Il2008[ ). Ee = 3.92 x 10^^ erg (Fig. [T]). The thin line shows the best fit with 
(3 = 1, whose TeV spectrum is not broad enough to explain the HESS observations. In 
the hard X-ray band, our model predicted emission spectrum is significantly harder and 
brighter, which gives better fit to the Suzaku observations and can be tested with future 
HXMT and NuSTAR observations. Our model also predicts a higher MeV fiux testable with 
future GLy4>ST observations. Figure [2] shows the dependence of the spectrum on p and (3. 

The distribution of relativistic electrons obtained above implies L = 2.34xlO^^[C2p/ {C2p+ 
l)]-i/2^o^i g^j^^ ^i^g electron density Ue = 2.00 x 10"=^(C2p + 1)5/4^^-3/2^-^1/2^-5/2 ^^.3^ 

where Y ~ 0.1 is the Helium abundance. For the steady-state solution to be applica- 
ble, Ci < 0.824(72^^ < 0.64. For C2p = Ci = 0.5, we have L = 4.04 x W'^Uq^ cm. 
He = 1.33 x 10~^Uo^^^ cm-3, Tac = USUq^ yrs, r^sc = ^2.3Uq^ yrs, and r] = 0.38A%Uo^^. 
The energy carried by the magnetic field Eb = = 2.77 x W^^Uq^ erg. The short 

coherent length of the magnetic field 2n/kd = 'jc'm'eC^ /qB = 1.02 x 10^^ cm is consisten t 



with the low level of polarization from most part of the radio image (ILazendic et al. II2004I ). 
If kd doesn't change significantly in the downstream, the distance that relativistic elec- 
trons diffuse through the t = 1600 yr lifetime of the remnant is about (t/T(,sc)^^'^L = 0.805 
pc. Most of them are therefore trapped near the SF. The total kinetic energy carried by 
the supersonic dissipative layer in the lab frame is ~ 27ri?^L(l -|- 2y)mpne(3f//4)^ = 
7.34 X IO^^Uq^^^ erg. The corresponding gas mass M = 2K/{?,U/Af = 1.29 x lO^^f/"^/^ 
g, thermal energy Ef ~ QuR'^Lmp'rif.il + 2Y)v'^ = 3.67 x lO^^f/g'^''^ erg, gas temperature 
ksTp = 2mp{l + 2Y)Et/3M{l — Y) = 26.4f/Q keV, where we have assumed that the temper- 
atures of protons and Helium ions are the same and Tg <C Tp. 

The acceleration efficiency of ~ is not sensitive to p and U and therefore doesn't 
change significantly through the evolution of the SNR, implying that a thermal energy of 
Ee/v ~ 1-02 X IO^^Uq^^"^ erg, which is 27.8 times higher than Et, be produced in accompany 
with the electron acceleration. The mass of the shocked thermal plasma is about E^M/ rjEt ~ 
3.58 X IO^^^/q g= 1.79?7o ^^^M©. These suggest that most of the relativistic electrons were 
accelerated more than Tac ~ 100 years before outside the current dissipative layer. If the 
shocked plasma is uniformly distributed within a shell of 4.5 pc (Aharonian et al. 2006), the 
electron density is 1.74 x lO-^f/,^^''^ cm"^, which is comparable to that inferred from X-ray 
observations (Cassam-Chenai et al. 2004). 
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However, X-ray observations of the north-west outer edge of the remnant suggests U < 
0.015c (Uchiyama et al. 2007). So the model predicted thermal X-ray emission likely exceeds 
the observed value. The north-west edge is interacting with dense molecular clouds. It is 
possible that the shock speed in a low density medium is still greater than the observed upper 
limit so that U > 0.015c. On the other hand, given the uncertainties in the coUisionless 
electron heating processes, the electron temperature may be well below 1.0 keV so that 
no significant thermal X-ray emission is expected in the XMM-Newton and Chandra X-ray 
bands. Coulomb col lisions will make fceTe ~ 0.198(ne/0.01cm-3)2/5(t/l600yr)2/5[/o/^ keV 



(iHughes et al. Il2000l ). Line emission dominates the cooling with a luminosity of ~ 10 erg 



s ^ and may be observed in the UV to soft X-ray band. Detailed modeling of the e volution 



of the SNR is needed to reach more quantitative results fiDecourchelle et al. Il2000l ) 



X-ray observations discover thin bright features varying on a timescale of about a year. 
The width of these features is about w ~ 4'.'0 ~ 6.0 x 10^^ cm, corresponding to a speed 
> 0.063c. Since the shock speed is much less than this value, the variabi lity has been 
attributed to cooling in very strong magnetic fields (jUchiyama et al. 1 120071 ) . According 
to our model, the decay may be due to the diffusive escape of relativistic electrons from 
bright regions. For electrons with a scattering mean free path of 2n/kd, the corresponding 
diffusive escape time from a region of size w is about 0.93 years. The X-rays are produced 
by electrons with energies greater than 'jcmeC^, whose scattering mean free path by fast 
mode waves should be longer than 2TT/kd implying a diffusion timescale of a few months. 
The formation of these bright features is then related to the inhomogeneity in the pre- 
shock medium, which is supported by the relatively higher level of radio polarization nearby 
(ILazendic et al. II2004J ). The fact that the brightening features are not located at the outer 
edge of the remnant can be attributed to the projection effect. 



We also apply the model to SNR RX J0852.0-4622. For = 1 kpc, R = 17 pc, we 



have B = 9.4 fiG, p 



2, 7cmec2 



2.94 TeV (7c = 5.76 x 10^), and the electron energy 



Ee = 1.01 X 10^*^ erg. These imply L 



3.73 X W^Uq^ cm, 



M 



141f/o ^ years, Tesc 
2.01 X 1032f/-^/2 



35.3f/o ^ years, rj 



0.404%[/o/^ 



7.76 X 10-3?7o ^''^ cm-3, 



K 



1.14 X W^U, 



-3/2 



erg. 



5.71 X W^Uq^''^ erg. 



26.4f/2 keV, and Eb 



= 4.53 x W^U^^ erg. These two remnants are very similar. 



4. Conclusions 



The steady-state distribution of relativistic electrons accelerated by fast mode waves in a 
downstream turbulent dissipative layer with a thickness L of weakly magnetized coUisionless 
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astrophysical shocks can be approximated as 

/ oc 7~^exp-(7/7c)^/^ 

where p = (9/4 + 2'^^3^c'^vy5^U^y/^ - 1/2, -fcmec'^ = qBL5^/^[{p + 0.5) - 2.25]i/2^/24c 
for Ci = C2p = 0.5, and the ratio of the dissipated energies going into nonthermal and 
thermal particles is ~ = (e^/^96/57r'^)^/^f^/f/^. For p > 2.09, the turbulence decay time 
is shorter than the particle acceleration time and one needs to consider the time dependent 
solution. The accelerated electrons in the STSNRs give excellent fits to the broadband 
spectra of SNRs RX J1713.7-3946 and J0852.0-4622. The model attributes the recently 
observed variable X-ray features to inhomogeneity in the upstream and the spatial diffusion 
of relativistic electrons near the cutoff energy and predicts harder hard X-ray spectrum, 
higher hard X-ray and 7-ray fluxes than leptonic models with a sharper cutoff in the electron 
distribution. Future high sensitivity hard X-ray and 7-ray observations will test the model. 
The major uncertainty of the model is related to the turbulence generation mechanism near 
the collisionless SF. In this paper, we assume that the turbulence is isotropic and has a 
characteristic eddy size of L. Little is known about the evolution of anisotropic turbulence. 
The latter is valid as far as the free energy dissipation through shocks proceeds over the 
scale L, instead of, instantaneously at a narrow region (<C L) near the SF as suggested in 
the diffusive shock models. 
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Fig. 1. — Model fit to the broadband spectrum of SNR RX J1713-3946. Tlie low and 
high energy spectral humps are produced through the synchrotron and IC scattering of 
the Galactic interstellar radia t ion, re spectively. The X-ray data are for the four models 
discussed in iTakahashi et al. I (|2008al) and rescaled to the ASCA flux level for the whole 



20031 ) . The thick solid line corresponds to the fiducial model 
7.2 X 10^ and B = 12.0/iG. The thin solid line with p = 2.0, 



remnant (lUchiyama et al. 
with (3 = 0.5, p = 1.85, 7c 
B = 12.0/iG, and 7c = 4.4 x 10^ is the best for an electron distribution with an exponential 
cutoff, i.e., P = 1.0. The fiducial model gives a better fit to the X-ray and TeV spectra. Insert 
is an enlargement of the TeV spectrum to demonstrate the improvement of the fitting with 
the fiducial model. The dashed line gives the GL^S'T integral sensitivity (1 yr 5a: McEnery 
et al. 2004). The dot-dashed and dotted lines correspond to the HXMT (lO^s 3a: Li et al. 
2006) and NuSTAR (lO^s 3(t: www.astro.caltech.edu/~avishay/zwickyl/ppts/Harrison.ppt) 
sensitivities, respectively. 
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Fig. 2. — Dependence of the emission spectrum on (3 and p. B = 12 G. 7c and are 
adjusted to fit the radio to X-ray fiuxes. The two dashed fines are for p = 1.85 and (3 = 0.25 
and 1.0. The thin sofid fines are for j3 = 0.5 and p = 1.7 and 2.0. 



